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Abstract 

The theory of string-like continuous curves and discrete chains have numerous important physical 
applications. Here we develop a general geometrical approach, to systematically derive Hamilto- 
nian energy functions for these objects. In the case of continuous curves, we demand that the 
energy function must be invariant under local frame rotations, and it should also transform covari- 
antly under reparametrizations of the curve. This leads us to consider energy functions that are 
constructed from the conserved quantities in the hierarchy of the integrable nonlinear Schrodinger 
equation (NLSE). We point out the existence of a Weyl transformation that we utilize to introduce 
a dual hierarchy to the standard NLSE hierarchy. We propose that the dual hierarchy is also in- 
tegrable, and we confirm this to the first non-trivial order. In the discrete case the requirement of 
reparametrization invariance is void. But the demand of invariance under local frame rotations pre- 
vails, and we utilize it to introduce a discrete variant of the Zakharov-Shabat recursion relation. We 
use this relation to derive frame independent quantities that we propose are the essentially unique 
and as such natural candidates for constructing energy functions for piecewise linear polygonal 
chains. We also investigate the discrete version of the Weyl duality transformation. We confirm 
that in the continuum limit the discrete energy functions go over to their continuum counterparts, 
including the perfect derivative contributions. 
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I. INTRODUCTION 



The theory of continuous curves in three dimensional space pQ is one of the pillars of 
differential geometry [2J. The study of string-like objects and their properties, both contin- 
uous and discrete, is similarly pivotal to several apparently disparate subfields of physics. 
Examples include polymers [B], fE], Kirchoff-type elastic rods [5]- [7], vortices in fluid dynam- 
ics, [8], [9], [10], turbulence [UJ, superconductors [12], [13], superfluids [Uj, plasmas [15] . 
metallic hydrogen [16], and confining [T7], cosmic [18] and fundamental [19] strings in high 
energy physics, and numerous other applications. Strings can display intricate geometry 
including knots [20], strings can proceed by leapfrogging [21], [10], and strings might even 
realize exotic exchange statistics in three dimensions [22]. Among the several general level 
theoretical contributions, we in addition draw attention to |23J and to |24J as particularly 
notable ones. 

In the present article we shall be mainly interested in deriving geometric energy functions 
for discretized curves, or more precisely piecewise linear polygonal chains in three space 
dimensions. The motivation comes from a recently proposed approach to the physics of 
proteins [25]- [29]. A protein is a biologically relevant example of a piecewise linear polygonal 
chain, with vertices identified as the central C a carbons of the amino acids. A biologically 
active protein is also particularly interesting from the point of view of physics of string-like 
structures in three space dimensions. It is an object that despite an inherent one dimensional 
character, has physical properties that causes it to collapse into a space filling construct. 
We propose that the study of space filling, collapsed string-like objects can have several 
important ramifications also in other areas of physics, from understanding turbulent flows 
in fluid dynamics to potentially new forms of matter in high energy physics. 

Our approach is based on the observation, first made by Hasimoto [30], [31], [H], [9] in the 
context of fluid dynamics, that the one dimensional nonlinear Schrodinger equation (NLSE) 
describes string-like objects such as vortex filaments. For this he introduced a change of 
variables that relates the wave function of the NLSE to the Frenet frame representation of 
a space curve. The NLSE is a widely studied example of an integrable model [32]- [M] both 
due to its stature as a universal theoretical construct, and due to the abundance of both 
theoretical [32], [33] and practical [3J] applications. Among the remarkable properties of 
the NLSE is the ability to support solitons as classical solutions. In the case of continuous 
curves, when the energy is computed by the canonical Hamiltonian of the NLSE and the 
ensuing configuration is related to the geometry of a curve by the Hasimoto transformation, 
the soliton describes the buckling of the curve. 

A properly discretized version of the NLSE preserves the integrability [35] , [36J . It also 
supports soliton solutions that can be utilized to model the buckling of piecewise linear 
polygonal chains, via a discretized version of the Hasimoto transformation. This has been 
applied in [27], [37] to propose that the loop structures that are the characteristic geometric 
features of folded proteins, could be modeled using a properly discretized version of the NLSE 
energy function. In the case of proteins, the loops are solitons that cause the string-like C a 
backbone to buckle and collapse into a space filling configuration. As a consequence, from 
a general point of view, chains that support solitons with geometric interpretation are truly 
curiosity arousing: Despite its inherently one dimensional character, a discrete string-like 
chain that supports solitons athat cause it to buckle, appears like a localized space filling, 
particle-like extended object. 

In the following, we shall first outline how energy functions of continuous curves are 
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derived from purely geometric considerations, by demanding both invariance under local 
frame rotations and covariance under curve reparametrizations. We show how the NLSE 
together with its integrable hierarchy emerges naturally, when we demand that the energy 
function of a curve can not depend on the choice of framing along the curve. We scrutinize the 
properties of the Zakharov-Shabat [38J recursion relation between the conserved quantities in 
the NLSE hierarchy, from the point of view of frame independence. We observe that there is 
the possibility to introduce a dual hierarchy, which is related to the original NLSE hierarchy 
by a Weyl scaling transformation. It turns out that the dual hierarchy has some quite 
attractive features. In particular, it embraces two of the conserved quantities of the NLSE 
hierarchy, the length and the helicity (Chern-Simons), that have no place in the standard 
NLSE recursion relations. We confirm by explicit computation that the leading order terms 
of the dual hierarchy are mutually in involution, suggesting that the dual hierarchy gives rise 
to a novel integrable set of conserved quantities. But at the moment we still lack a general 
proof that the dual hierarchy is fully integrable i.e. that all the infinitely many quantities 
that form this hierarchy are indeed in involution. 

We then proceed to the discrete case, of potentially space-filling chains. We first review 
the discrete generalization of the Frenet equations [39]. We show by an explicit compu- 
tation, that the two known integrable discretizations introduced in [35] and [36] coincide 
[30] , by devising an explicit change of variables that maps these models to each other. We 
then proceed to a discretization of the Zakharov-Shabat recursion relation. For this we 
adopt a geometric principle, that the discrete version of the recursion relation should be 
formulated by demanding invariance under local frame rotation. We introduce the ensuing 
discretizations of the conserved quantities, both in the case of the standard NLSE hierarchy 
and in the case of its Weyl dual hierarchy. We confirm that in the continuum limit where 
the distance between the vertices in the discrete case vanishes, the discretized versions of 
the conserved quantities indeed go smoothly over to their continuum versions, including the 
perfect derivative contributions which are necessary for the consistency of the higher order 
Zakharov-Shabat recursion relations. We conclude with a discussion on the discrete energy 
functions that can be utilized for example in numerical studies, to describe the physical 
properties of piecewise linear polygonal chains, their buckling, and the ensuing space-filling 
capacity. 



II. FRAMING AND REPARAMETRIZING CONTINUOUS CURVES 

We consider a general, filamental, space curve 7 in M. 3 that can be described by a class 
C 3 differentiable map x(z). Here z G [0, L] is a generic parametrization and L is the length 
of the curve, 

L L 

H i = L = J dz y/x z ■ x z = J dz\\-^\\ (1) 


We note that this is essentially the standard time independent Nambu-Goto action of a 
string. For simplicity we shall assume that there are no inflection points along 7, so that 
||x z || 7^ 0. But a generalization to include inflection points is straightforward |39j . 

In the absence of inflection points, the curve 7 can be globally framed as follows. The 
unit length tangent vector is 

* = m^ x - ( 2 ) 
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It is orthogonal to the unit length bi-normal vector 



| |X Z X X Z2 1 1 

The unit length normal vector is 

n = b x t 

The three vectors (n, b, t) form the right-handed orthonormal Frenet frame, at each point 
of the curve 7. 

The Frenet equation relates the frames at different points along 7. Explicitely, the Frenet 
equation is PQ, [2] 

'n\ /Or -k\ /n N 

V9 -r b I (3) 




Here 



v 7 ^ = 3" = V*z(z) ■ xjz) = \\x g \\ = ^fgz~ z (4) 



dz 

Since 7 is one dimensional, the metric on 7 has only one component that we denote g zz . 
This determines the local scale of curve parametrization vis-a-vis the metric in M. 3 . Further, 

< Z ) = || x 113 ( 5 ) 
1 1 z 1 1 

is the (Frenet) curvature of 7 on the osculating plane that is spanned by t and n, and 

z \ \X Z x yi zz ) ■ x 222 
T(Z) = llx xx IP (6) 

I \ A z ^ A Z2 1 1 

is the torsion of 7. 

If the metric, curvature and torsion are known we can construct the frames by solving 
([3]), and then proceed to construct the curve by solving ([2]). The solution is unique, up to 
rigid translations and rotations of the curve. In this manner the Frenet equation proposes 
us to construct energy functions for curves in IR 3 , in terms of the curvature and the torsion 
and the way how the curve is parametrized. 

In the following we shall denote by s G [0, L\ the arc-length parameter, while z denotes 
generic parametrization. The arc-length parameter s measures the length along 7 in terms 
of the distance scale of the three dimensional ambient space. The change of variables from 
a generic parameter z to the arc-length parameter s is 



s(z 



\x z {z')\\dz' 







Accordingly, we consider the effects of infinitesimal local diffeomorphisms along 7, obtained 
by deforming s as follows 

s -)• z = s + e(s) (7) 
Here e(s) is an arbitrary infinitesimally small function such that 

e(0) = e(L) = = e s (0) = e s (L) 
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The Lie algebra of diffeomorphisms ([7]) of a line segment in IR 1 is the classical Virasoro 
(Witt) algebra. This proposes us to define a function f(s) on 7 to have a weight h akin the 
conformal weight, if f(s) transforms according to 

6f{s) = -(e± + he.'\ f(s) (8) 

under the infinitesimal diffeomorphism ([7]). If / has weight h 1 and g has weight h 2 then the 
product fg has weight hi + h 2 . We also note the finite version of (JsJ) , 

r// x 



Since the three dimensional geometric shape of the curve 7 in IR 3 does not depend on 
the way how it has been parametrized, the embedding x(z) transforms as a scalar i.e. it 
has weight h = under reparametrizations. Similarly, the curvature ^ and torsion ^ are 
scalars under reparametrizations. Infinitesimally, 

8k(s) = -e(s)-^ = -eK s 
as 

5r(s) = -e(s)^ = -er s 
as 

Generic functions can be composed by taking derivatives, multiplying and summing up 
functions that have definite weights. 

If f(s) is a function on 7 that has a definite weight, in general its derivative does not have 
a definite weight. For the derivative to acquire a definite weight, we need to extend it into 
a covariant derivative along 7. Given the metric g zz along 7, we can deduce the covariant 
derivative by considering its action on a function / with weight h. For this we demand that 
under the infinitesimal transformation (|7|) 

s +r w) /w_ (s +rw ) /w 

-(e4- + (h+l)e,) (±+r(,))f(s) 



ds J \ds 

Here T is the connection on 7. Explicitely, we can choose 

r{z)=hd z hiy/g = hd z ln\\x z (z)\\ (9) 

where we have used Q. More generally, we may also choose 

T(z) -)• T(z)+A(z) (10) 

Here A(z) is a quantity with weight h — 1. 

A general frame is related to the Frenet frame by a local SO(2) rotation, that sends 

n \ (^ e A = f cos? ?(» -sin?7(2;)\ /n\ 
by \ e 2/ V sin 77(2) cosrj(z) J I by 
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The ensuing rotated version of the Frenet equation is 








r -Vz 




V9 I - r + Vz 

k cos rj k sin 7] 



-kcostA I ei 
-Ksin?7 I e 2 

o / V t 



We can write this as 



d_ 

dz 



ei 
e 2 



T r 

( -T r -K„ 






(12) 



(13) 



Here k 9 is the geodesic curvature, and K n is the normal curvature, and r r is the relative torsion 
of the curve on a two dimensional surface 5 G M 3 ; we obtain the surface by deforming the 
osculating plane, around the point of contact with the curve, and different values of T] can 
also be interpreted in terms of different choices of the surface S. Notice that while n g and 



n n in Q13|) both depend on the surface, that is they are both frame dependent, the modulus 

(14) 



P 



is a frame independent characteristic of the curve 7. In particular, p is a scalar under 
reparametrizations. 

We introduce the combination 



ip{z) = K,(z)exp J i J Ty/gdz' 


Alternatively, we may introduce the more symmetric 



(15) 



M :K'X|> I - jTy/gdz' +- I r^fgdz' 



The decomposition ((15j) is essentially the Hasimoto variable of fluid dynamics 

is invariant under 

t(s[z\) -> t(s[z]) 



has weight h = 0, and it is invariant under the frame rotation (12): When 

drj 



ds 



r{z) 



1 drj 



y/gdz 



and 

we get 



k ->■ e vq K 



ijj(z) -> k e tv ^ ■ exp I i I </g r — ir](z) ) =ip(z) 



(16) 
. It 

(17) 

(18) 
(19) 



We also note that there are the following two natural, frame independent realizations of the 

(20) 



quantity A with weight h — 1 in (10) 



■A rsj j z (z) = d z ]nijj 



and 



A ~ j z (z) = d z lnip 
They both appear in our construction. 



(21) 
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III. NLSE HIERARCHY OF CURVES 



We are interested in frame independent energy functions of curves, that have definite 
and identifiable transformation properties under curve reparametrizations. We first observe 
that there are several frame independent energy functions, that have a natural geometric 
interpretation. The length if_i = L in is an example. Additional familiar examples are 
the total torsion (total helicity) 



H_ 2 = / r(s)ds (22) 



and the total squared curvature 



L L 



Hi = \J K 2 (s)ds= l -J \t s \ 2 ds (23) 



The latter determines the Worm Like Chain (Kratky-Porod) model that has been extensively 



applied e.g. to study elastic properties of DNA [H]; for applications of (23) to random 
surfaces, see [21]. The Bernoulli elastic curve is defined by the following linear combination 
of ([23} and Q, 



Hi + tf-i = \ /{k 2 (s) + A}J.s 



o 



Here A is a Lagrange multiplier that enforces the length constraint. Finally, we also mention 
the following combination of curvature and torsion 



L 



H K = l -j{a k 2 (s) + (3 r 2 (s)}ds (24) 
o 

This defines the Kirchhoff elastic energy of a filament, that has also been studied extensively 
in the literature [5]- [7]. Here we have assumed the arc- length parametrization; recall that 
both r and k are scalars under reparametrizations. 

In the context of three dimensional fluid dynamics, the length is commonly employed 
as the Hamiltonian of a vortex filament [5]-[TT]. The dynamics of the filament is given by 
the localized induction approximation 

dx 

— = x s x x ss = Kb (25) 



where the time evolution follows from the Rasetti-Regge bracket [23] . with symplectic one- 
form 

ds dt x ■ x s x xt 



Similarly, the energy (|23|) has the form of standard Hamiltonian of Heisenberg spin chain, 

ven by 

t x t„ (26) 



with t(s) the spin vector. The dynamics is given by 

dt 

dt 
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and now the Poisson bracket is 



{t a (s),t»(s')} 



abc±c 



t c (s)6(s 



(27) 



Hasimoto [30], [31 J utilized the variable (15) to establish an equivalence between (25) and 
(26). For this he showed that both are equivalent to the equation of motion of the non-linear 



Schrodinger (NLSE) Hamiltonian 



ds + A(^) 2 } 



when the Poisson bracket is 



{ip(s),i/;(s , )} s = i5(i 



(28) 



(29) 



Here we have introduced a generic coupling parameters A, it can be scaled to A = 1 by 
suitable redefinitions. 



In the following we adopt (29), its representations in terms of the Hasimoto decompo- 



sition, and its discrete version as the standard bracket that we denote with the subscript 
s. 



The NLSE Hamiltonian system (28), (29) can be elegantly presented in a form which is 



manifestly covariant under reparametrizations of the curve: When we introduce an arbitrary 
parametrization s — > z, the NLSE Hamiltonian becomes 



J-gdz {g zz ^ z + X(^) 2 } 



(30) 



where we have used the definition Q of the induced metric on 7. The Poisson bracket is 
similarly reparametrization covariant, 



{iKz),$W}. = —6{z-z>) 



(31) 



The NLSE equation is known to be integrable on the finite segment [0, L], when ap- 
propriate boundary conditions are imposed. In fact, the equations (25), (pBl constitute a 



bi-hamiltonian pair that determines the NLSE hierarchy. Consequently there is an infinite 
set of conserved charges in involution with respect to the ensuing brackets, 



{H n ,H m } = n,meZ" 
The charges are integrals of densities 



(32) 



H,„ 



y/gdzUniz) 



that are functional of if), ^ and their derivatives. 

In general, we can expect that the densities in an integrable hierarchy are determined 
only up to an additive derivative contribution, 



n n ~ H n + dk 



(33) 
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where A(z) is some functional of ip and ip and their derivatives, and subject to an integrable 
boundary condition such as A(0) = A(L). Recall that in one dimension, any function can be 
represented as a derivative of its own integral. Consequently, we limit our attention in (33) 



to such functional A(z) that are local in the fields ip(z) and ip{z) and their derivatives. In 
that case, we call dA(z) a perfect derivative. In an integrable hierarchy the density T-L n {z) of a 
conserved charge is in general an equivalence class, where two representatives are equivalent 
to each other when they share the same Poisson bracket relations and deviate from each 
other only by a perfect derivative. This equivalence relation defines a structure akin de 
Rham cohomology. A mutually consistent choice of the perfect derivative contributions in 
a hierarchy of conserved charges in involution, can be viewed as a "choice of gauge". 

Our aim is to utilize the NLSE hierarchy to introduce frame independent energy func- 
tions, in a manner which is manifestly reparametrization covariant. Occasionally, we find it 



convenient to introduce j z (z) in (20) and 



p(z) = if>(z)ij)(z) 

as the canonical variables, in lieu of the Hasimoto variables ip(s) and ip(s). Note that 
alternatively the complex conjugate j z (z) in (21) could be used instead of j z (z), and even 
a linear combination of the two could be utilized. In terms of the variables j z and p, the 
NLSE Hamiltonian has the following manifestly covariant representation 



yfgdz { g zz j z p z - pg zz 3zj z + Xp 2 } 



(34) 



and the standard bracket (\31\j is the inverse symplectic two-form 





n-l(z,z') 



{jz (z) , jz(z') } s {j z (z) , p(z') }, 
{p(z),jz(z')}s {p{z),p(z')} s 



d 1 

dz ^/g 



1 d 

y/g dz 





5(z - z') (35) 



We also record the ensuing NLSE equations of motion, in terms of the arc-length parameter: 

d 



It J = [H ^ j{s)}s = l Ts 



ds 



j-f-2Xp 



and 



, d 
"ds 



(36) 



(37) 



We remark that both (36) and (37) are perfect derivatives. Here and in the sequel we denote 

3 = 3 s 

the variable j z (z), in arc-length parametrization. 



Instead of (25) and (26), it is often more convenient to base the bi-Hamiltonian structure 



on (34), (35), and with the conserved momentum 



P = H 9 



dsj p 



Vgdz g zz j z p 



(38) 



{# 3 ,# 2 }s = 
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as the second Hamiltonian, together with if 3. In the arc- length parametrization, the sym- 
plectic structure corresponding to H 2 is 



{p(s),j(s')} {p(s),p(s')} 



2A# 

as 



ds^dsJ > 



d -M -(£ P + p d 



5( S - S) 



(39) 



• ds 



ds ■ 



It is straightforward to confirm that with (38), (39) we obtain (36), (37) as the equations of 
motion. The infinite number of conserved densities W n (s), n = 1, ... of the NLSE hierarchy 
can be iteratively constructed as follows, 



/ sn n+1 \ 

' Sj(s) 1 



SHn 

Sj(s) 

SHn 

S P (s) 



(40) 



and the construction is conventionally started from the conserved charge ( 23 ) 

ds k 2 



#1 



dsK\ 



dsipip 



(41) 



We note that the length ([I]) and the total torsion (22) are also conserved charges in the 
NLSE hierarchy, even though they do not have a natural place in the standard hierarchy 
relations; the square of the torsion that appears in (24) is not the density of a conserved 
charge in the NLSE hierarchy. As a consequence, both ([I]) and (22 ) are commonly interpreted 
as additional conserved charges, with negative order H_i and H_ 2 respectively (see e.g. [5J, 
[S]), and included ad hoc. We may also choose 

#o = 

Finally, we point out the following relation between the NLSE hierarchy and the classical 
Virasoro (Witt) algebra of diffeomorphisms on the segment [0, L] eR 1 . If we define 



Pis) 



-. 00 



e 2 ™ n z 



(42) 



and 



j 00 

J2 5Z^ r ' 



and substitute in (39) the coefficients obey the Poisson brackets 

cL 



2ii(n — m)C 
—1 r Km^J rt j rm 



{^ni <^Jm\ s 
\\Jni JmJ s 

{L,£ n } s = {L,^T n } 



i 4Tr 2 n 2 5 n _ m L 



47rAn<L _ m L 



(43) 



(44) 

(45) 
(46) 
(47) 
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where (44) is the classical Virasoro (Witt) algebra and J n extends it with a current algebra; 
the central charge c is classically absent but in general it is non-vanishing, in the quantum 
theory. Note that we have here interpreted the length as an additional element of the 
algebra; it is after all a conserved charge in the NLSE hierarchy. We have also renormalized 



the coupling A in (28) as follows 



LA = A 



Moreover, 



Hi 



dsipip 



dsp oc C c 



where we use (42). Since ^3 in (35) is field independent and can be presented in terms 



of the current algebra generators (44)-(46) we conclude that all the classically conserved 
charges H n are polynomials of the classical Virasoro-current algebra generators C n and J n , 
and length L. 



IV. REPARAMETRIZATION COVARIANCE OF RECURSION RELATIONS 



Besides (40), the NLSE hierarchy of conserved densities / H n {s) with n G Z + can be 
constructed using the following recursive relation (38], [32] 



u n+ i(s) = -i 



du n 
ds 



n-l 
3=1 



(48) 



We assume that the arc- length parametrization is used, in the case of curves in IR 3 . We 



also note that alternatively, the recursion could be based on the complex conjugate of (48). 



When the construction in (48) is started with 



u {s) 
^i(s) 





|V>(s) 



(49) 



we reproduce the densities H n (s) of all the higher order conserved charges as follows 
[33], 

■H n {s)=ipu n neZ + (50) 
For a generic parametrization s — > z of the curve, we need to covariantize the recursive 



relation (48), with due care. For this we remind that the Hasomoto variable ip(z) is a scalar 
i.e. has weight h = under reparametrizations. Thus we obtain the reparametrization 



covariant expression of the charge Hi in (23) by interpreting the ensuing density coi as a 
scalar. In a generic parametrization, 



Hi 



^fgdz ui 
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Next, we aim to compute £^2(2:) from the proper reparametrization covariant version of (48). 
Since u\ is a scalar, we conclude that we may set 



u 2 = -t—cui 

dz 



1. In (38) we have displayed the corresponding conserve charge 



so that u 2 has weight h 
if 2, in a covariant form. 

We proceed to H 3 . Since u\ has weight h — 1, for a covariant charge we need to interpret 
the derivative term in (48) accordingly, as the weight h = 1 covariant derivative using the 
connection d9l). The ensuing covariant version of (48) then reproduces (30): 



u 3 



^31 + w 32 



d 

dz 



+ T)g zz uj 2 + XiI)ujI 



=> H 3 = J^dzi)u 3 = i J^/gdz {g zz i) z i) z + A(^) 2 } 

Note that in u 3 (z), the first term u 3 \ has weight h = 2 while the second term CJ32 has weight 
h = 0. Consequently, when we proceed to if 4 and onwards, and compute the derivative 
of Co>3 etc. in the recursive relation (48), the derivative operator must be replaced by the 



appropriate covariant derivative: In the case of f/4 we need to introduce a covariant derivative 
with weight h = 2 when the derivative is acting on the first term CJ31, and a covariant 
derivative with h = when acting on the second term in u 32 , 



d 
ds 



-i(d z + 2T)u 31 - id z u 32 



In this manner, by successively interpreting the derivative in (48) as the appropriate co- 



variant derivative, we generate the manifestly reparametrization covariant versions of the 
densities u; n+ i(z), order by order. These densities are linear combinations of terms with 
different weights, and the highest weight term has weight h = n. 



V. HIERARCHY WITH POSITIVE ORDER 



We first remind that the densities 7f n (s) are not necessarily uniquely determined. In 
particular, there may be the latitude (33) to add a perfect derivative to the density. In 



the case of an integrable model, this latitude can be restricted by imposition of proper 
(integrable) boundary conditions at the end points of the segment [0,L]. Here, our goal is 
to utilize the NLSE hierarchy to find frame independent and reparametrization covariant 
energy functions for curves. For this we note, that the end points of our curves can move 
freely in M. 3 . As a consequence the boundary conditions that we need to impose at the end 
points z = and z = L, are to be open. In particular, we may deform the energy density of 
the curve by an appropriate perfect derivative term, as long as this term does not interfere 
with frame independence and reparametrization covariance. 

From the point of view of algebraic relations in integrable models the canonical pair 
of complex variables (tp, ip) is adequate. But we are interested in describing curves. For 
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this, more work is necessary. We need, in addition of solving the integrable equations, 
to construct a solution to the Frenet equation ([3]), (13), followed by an integration of the 
curve reconstruction equation (|2]). For this reason we proceed to investigate the recursion 
relation (48) explicitly, in terms of the (Frenet frame) curvature and torsion. For clarity, 
in the sequel we use the arc-length parametrization, a reparametrization covariantization is 
straightforward following Section IV, and will be performed when conceptually desirable. 



We first observe that the relation (50) between the density T-i n and the quantity u n admits 



the following internal symmetry: The density remains intact, when we send 

u n {s) ->■ e ia(s) w n (s) 

(51) 

where a(s) is some function. Note that a(s) can also depend on the canonical variables. In 
particular, (51) does not need to be a canonical transformation. The transformation (51) 



converts the recursion relation (48) into 



n-1 



d da\ ■sr^ 

% ~ds + ~ds ) Un+ ^2-^ ^j^n-j 
j=i 



We recall the Hasimoto variable ( 15 ) in the arc-length parametrization 



(52) 



ip(s) = K,(s)exp I i / r ds' 



with (k, t) the Frenet curvature and torsion, and choose 



a(s) 



fs) ds' 



This gives us 



Wn+l(Sj 



- I 1-. T 

ds 



\ n—1 
» J U n + \4> UJjUJn-j 
' .7 = 1 



(53) 



which is manifestly covariant under frame rotations. 

We record the following five terms, for later reference; in computing these, we employ the 



Frenet frame a.k.a. unitary gauge relation (48), in which the curvature k(s) is a real valued 
quantity. 

uo = 



4« 

\tk-\k' 



-\k" + \t 2 k+ A K 3_i( r / K + 2rK /) 



A ( 4/tr 3 - Akt" - 12tk" - 12rV + 3Xk 3 t ' 



(54) 
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'4/c" 



Utt'k - 12k't 2 - 5\k 2 k') 
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We multiply these with R to obtain the ensuing versions of the integrable densities r H. n \ we 
remind in the Frenet frames the curvature is real so that R = k. We list the following four 
for later reference: 

Hi = y 

H 3 = \Wf + \t 2 k 2 + A K 4 _ id {TK 2 _ iKKl) 

^ = \kk'" (55) 
+ ^ ( 8kk"t + Ak 2 t 3 - A(k') 2 t + 3Ak 4 t ) 
+ ~{ -2Ak 2 t 2 + 5Ak 4 + m [ k 2 t' + kk't]) 



We identify in Hi the mass (number) density of NLSE. In I-L2 we identify the canonical 
momentum density, in addition of a perfect derivative term. In we have the NLSE 
Hamiltonian density, in addition of a perfect derivative term. Finally, H,^ is the complex 
modified KdV density together with a perfect derivative term. Note that the form of the 
perfect derivative terms is fully dictated by the recursion relation. Even though the presence 
of a perfect derivative does not contribute to the charge at the given order, it does influence 
the functional form of higher order conserved charges. 

We can also present the recursion relation directly in terms of the integrable densities. 
The result is 

H =0 



1 ^2 

(56) 



V =i + (ir- M) = £ - £ln^ = £-j(s) 



It is notable that the h = 1 quantity (21) appears. 



VI. HIERARCHY WITH NEGATIVE ORDER 



We observe that the derivative operator V that appears in (56) and its conjugate 



are related by the Weyl transformation, 

V = e - 9 Ve 9 (57) 

when we choose 

e e = H 2 (58) 
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This symmetry between the two quantities pOl) and (21 ) is most natural, and proposes us to 



introduce the following extension (Weyl dual) of (48), for all negative integers: We simply 
map 



This gives 



with 



W_( n +i) 



— I 



du- 



ds 



W-1 



1 

5 



We define 

"H_ n = 4ju_ n 

to arrive at the following recursion relations of densities 

n = o 



V 



iVH- n + A|^| 2 H-jH-^j) 



(59) 



ds 



where 



(60) 



We present the following three negative order densities, for reference and scrutiny, using the 
Frenet frame: 

U -^ S ) = -^ ln V> = T-i|; In K 



H-As) 



-^ln^-(|ln^) 2 + A|^| 2 
[r - i(ln k)'] 2 + 1 + Ak 2 - i [r + i(ln « 



(61) 



Remarkably, the two conserved charges ([!]) and (22) of the NLSE hierarchy now appear natu- 
rally, as quantities in the negative order recursive relation. We remind that in a conventional 



approach to NLSE, these two charges do not arise from the recursion relations (48). Instead 



they are introduced in an ad hoc manner [5J, [5] as additional charges. Furthermore, the 
density "H_3 can be given a familiar interpretation in terms of quantities that relate to gauge 
theories. This density is like the gauge invariant Proca mass for r, when we interpreted r as 
a gauge field and In k as the ensuing Stueckelberg field; the combination is gauge invariant 



under the analog U(l) gauge transformation (17), (18). See also next Section. 

Unlike Q and ( [22] ) , the integral of "H_3 is not a conserved quantity in the NLSE hierarchy. 
For example, we find, using the standard bracket (29) 



{H 3 ,H_ 3 } s = 2i ds 



2A|^| 2 (ln^)" 



^0 
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However, the iterative relation (59) is suggestive that the negative order densities "H_ 



determine an independent hierarchy of conserved charges that are in involution amongst 



themselves, and possibly including Hi and "H 2 , m terms of the Poisson bracket (29). 







for n, m > 



(62) 



We have confirmed this to the densities given in (61), for example we find 

f d 

{#_ 3 ,#-4} s = 3A 2 ds—iW) 



ds 







At the moment we lack a simple proof that all the negative order charges are in involution, 



the integrability (62) is a conjecture; the fact that the first four charges are in involution 



strongly suggests that the conjecture could be correct. 

In Figure 1 we summarize schematically the positive and negative order hierarchies, and 
the correspondence for the quantities that have a familiar interpretation, together with their 
Weyl duality relations. 

H ° \ 

V \ % 

Weyl 

Mass -4 ► Length 



Z-S 



I 



I Weyl 
Momentum ^ y Helicity 

y Weyl 
AHM/NLSE < >• Proca 



I 



Z-S 



Z-S 



Z-S 



I 



Z-S 



FIG. 1: Schematic relations between those conserved densities that have a natural interpretation 



in terms of gauge field theory, in the positive order (55) and negative order (59), (61) hierarchies. 



AHM/NLSE is for Abelian Higgs Model a.k.a. nonlinear Schrodinger equation, Weyl denotes 
relation by the Weyl duality between the hierarchies, and Z-S stands for Zakharov-Shabat recursive 
relation. Both hierarchies amalgamate at Hq = 0. 



VII. FRAME INDEPENDENCE AND GAUGE STRUCTURE 



We employ the recursion relations (56), (59) to specify Hamiltonian energy functions, 
that describe curves and their dynamics in KF. The explicit recursion relations that we 
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have presented utilize a fixed variable s G [0, L\ that we have chosen to identify with the 
arc-length parameter of the curve. However, we have also explained how reparametrization 
covariant generalizations can be introduced, order by order in the recursion. 

Notwithstanding, in order to describe a curve in M 3 we need in addition to solve the 
Frenet equation ([3]), (13). The Frenet equation involves the curvature and torsion, and as 
such it is frame dependent. We proceed to complement reparametrization covariance with 
manifest frame independence. 



According to (17), (18) the curvature and torsion (k, r) transform exactly like an Abelian 
Higgs multiplet (</>, A) under frame rotations, provided we identify the Frenet frame curvature 
k as the modulus of a complex (Higgs) scalar field <f), and the torsion r as the one dimensional 
U(l) gauge field A. In the arc-length parametrization, the action of frame rotation can thus 
be summarized as follows: 



Ms 



A(s) ~ t(s) -> t(s)-t) s (s) 

Consequently frame independent Hamiltonians for curves can be constructed in parallel with 
U(l) gauge invariant Hamiltonians of the Higgs multiplet (0, A). An example of the latter 
is the standard Abelian Higgs model Hamiltonian 



H 3 = jds ||(A + ^ )0 |2 + A | 0r 
If we identify the Hasimoto variable as follows, 

s 

ip(s) = cj)(s) exp{i / A} 



(63) 



(64) 



we recover the NLSE Hamiltonian (28). On the other hand, if in (63) we introduce the 
change of variables 



k — <f> — y a e 



to 



r ~ A^ J = A-|A ln0+ ^ ln( 



(65) 



r + f g arg k 



the NLSE Hamiltonian becomes 

#3 = 



(66) 



Like ip in (64), both J and o are manifestly frame (gauge) independent variables; J is 



commonly called the supercurrent variable, in the context of Abelian Higgs model; a is the 
Higgs condensate. 

The entire NLSE hierarchy can be similarly presented in terms of the frame independent 



variables (65). For example, for H 2 we get from (38) (note: a = yfp) 



Ho 



dsa J 



(67) 
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In the Abelian Higgs model, this is the canonical momentum; for H\ we have in ( 56 ) 



ds a* 



(68) 



which is the Higgs mass term, in the Abelian Higgs model. 



For the negative order density "H_2 we get from (61) 

. d 



The charge is 



t — i — In k 

ds 



H 



d d , 

r + — arg k — i— In a 
ds ds 



ds J 



(69) 



We identify (69) as a one dimensional version of helicity, or Abelian Chern-Simons term, for 



the gauge invariant variable J. For the density of H_ 3 we get similarly 

2 



H 



din k\ . d 
ds J ds 



d\n k 



T — l- 



ds 



+ 1 



-r + i—j + i 

ds 



(70) 



The conserved charge is 



#_3 = J dsJ 2 + L (71) 
where L is the length of the curve. In the Abelian Higgs model, we identify the first term 



in (70), (71) as the one dimensional version of the gauge invariant i.e. frame independent 



Proca mass, for the variable J. 



VIII. DISCRETE CHAINS 

In the discrete case of a piecewise linear space polygon i.e. chain, the reparametrization 
i.e. diffeomorphism covariance is not of concern. However, the requirement of invariance 
under frame rotations persists. Consequently, we start by revisiting the construction of 
frames for a piecewise linear polygonal chain in IR 3 |39j. We take ( % — 1,...,N) to be 
the coordinate sites of the vertices, and at each vertex site % we introduce the unit tangent 
vector 

U = , rm ~ ri , (72) 

| r i+l — r i\ 

We frame the chain by introducing at each vertex i the unit binormal vector 

bi = ti ~ 1 X U (73) 

|tj_x x tj| 

and the unit normal vector 

n, = b l xt i (74) 
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The orthogonal triplet (rij,bj,tj) defines the discrete Frenet frame at the vertex r» of the 
chain. The bond angles along the chain are 



«i = = arccos (t i+1 ■ t*) (75) 

and the torsion angles are 

n = r i+1) i = sign{bi_i x bj ■ tj ■ arccos (b m • b*) (76) 

If these angles are known, we can use the discrete Frenet equation 

cos k cos r cos k sin r — sin /t N 

- sin r cos r J I h { I (77) 









b, 


hi 


H 









sin k cos r sin re sin r cos re 




to construct the frame at vertex % + % from the frame at vertex %. Once we have the frames 
we get the entire chain, up to global rotations and translations, using 

fe-i 

r fe = ^ \r i+1 - r f | ■ ti (78) 

i=0 

This amounts to an integration of (|2]), in the discrete case. For simplicity, in the sequel we 
assume that the bond lengths i.e. distances between the neighboring vertices are constant 

|r i+ i = a (79) 

With no loss of generality we can set ro = 0, and choose to so that it points along the 
positive z-axis, so 

r - ! = az 



In parallel with the continuum curves, (77) proposes us to construct energy functions for 
discrete chains in terms of the bond and torsion angles. 

We note that (78) does not involve the vectors rij and bj. Thus, as in the case of 
continuous curves, we may introduce an arbitrary SO (2) rotation between them, without 
affecting the chain. At each vertex site % this frame rotation acts as follows, 

cos Aj sin A, N 

-sinAi cosAi | ( b | (80) 
1 

Here the Aj are the arbitrarily chosen local rotation angles. On the bond angles, the effect 
of frame rotations can be presented as follows, 

KiT 2 -> ^(cosAj+i^ + sinA.+iT 1 ) (81) 




This is the real SO(2) version of the complex U(l) rotation (18). Thus we may consider k% 
as a complex variable so that its real part coincides with the T^^component of (81), and the 
imaginary part is the T 1 component of (81). The frame transformation property is 

Ki ->■ e iAi+1 Ki (82) 
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The bond angle (75) is the modulus of this variable. On the torsion angles, the frame 



rotation has the following effect [39] 



Ti -> n 



Ai+i + A< 



(83) 



^4 priori, the fundamental range of the bond angle is [0, tt]. For the torsion angle the 
range is 7j G [— 7r, tt). Consequently we may identify Tj) with the canonical latitude and 
longitude angles of a two-sphere S> 2 . However, in the sequel we find it useful to extend the 
range of Ki into [— 7r, tt] mod(2ir), but with no change in the range of Tj. We compensate for 
this two-fold covering of § 2 by introducing the following Z 2 symmetry 



- TT 



for all j > 



(84) 



This is a special case of ( 80 ) , with 



A, 
A, 



7T 





for i > i + 1 
for Z < i + 1 



IX. THE DISCRETIZED NONLINEAR SCHRODINGER EQUATION 



Two a priori different integrable discrete versions of the continuum NLSE have been 
presented, the one in [35] and the one in [36J. We review these models, before we proceed 
to derive a discretization of (53), (56), using the frame independence of a discrete chain as 



our guiding principle. Following [32], we denote the two discrete models by LNSi and LNS2, 
respectively. It can be shown, that LNSi is equivalent to the lattice Heisenberg spin chain 
model (LHM) [32] ■ In [10] it has also been shown, that LNSi and LNS 2 are gauge equivalent 
the sense, that their Lax pair representations can be gauge transformed to each other. But a 
proof based on gauge equivalence of the Lax pair does not give an explicit relation between 
the variables. In particular, it does enable us to reveal, how to describe discrete curves, 
transferably in terms of the two versions of discretized non-linear Schrodinger equation. For 
this reason, we show how to establish the equivalence of LHM and LNS2 in terms of an 
explicit change of variables. 

We start by describing the lattice Heisenberg spin chain model i.e. LHM. With a three 
component unit vector at site % = 1, ...,N, the Hamiltonian is (for notational simplicity we 
set bond length a = 2 in (79) and take the classical spin parameter s in [32] to have value 

8=1) 



N-l 



H 



LHM 



- ^ ln(l + Sj • Sj+i) 



(85) 



and the Poisson bracket is the discrete version of (27) 



The equation of motion for Sj is the Landau-Lifschitz equation 



d,Si 
~dt 



{Hlhm, Sj} 



dH 



-Si x 



LHM 



(86) 



(87) 
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It reproduces (26) in the continuum limit a — > 0. 

We may use (85) to introduce an energy function for discrete polygonal chains. For this 
we simply identify the spin variable Sj with the tangent vector tj, defined in (72). Note that 
since only Sj ~ tj appears in (85), (86) the energy function (85) determines chain dynamics 
which is manifestly invariant under the local frame rotations (80). In terms of the discrete 
Frenet frame, the equation of motion (87) becomes 

1 



dtj 
dt 



} ^tan b i+1 - tan — b, 



(88) 



so that for our polygonal chain we obtain, from (72), the discrete version of the localized 

dri 



induction approximation (25), 



dt 



tan — b; 
2 



(89) 



Since the Landau-Lifschitz equation is integrable there are conserved quantities, equal in 
number with the canonical coordinates (or momenta). In particular, we note that the 
equation of motion (88) preserves the closure condition of the chain 



d_ 
dt 



N 



5> 



i=l 



This is the discrete version of the continuum statement that H_i in Q, ( [59] ) is a conserved 
quantity. 

In terms of the Frenet frame bond (75) and torsion (76) angles, with (Kj +lji , r i+lj j) = 
(/«i,Ti) the Poisson bracket (86) is 



{/%_!, Tj} 
{K i+1 , Ti} 



sm T i+1 

COSTi-i CSC 



COt + COS Ti COt Hi 



— COt y ~ cos T i c °t 
COS Ti+l CSC Hi 

{Tj_i,Tj} = CSC Hi— i (sin Tj COt Ki -\- sin Tj_l COt Kj— 2) 

{ri-i,T i+1 } = sinrjCSCKj_iCSC/tj 



(90) 



The Hamiltonian is 



H = -^ln(l+tj-t i+1 ) = -2^1n 



COS 



(91) 



and the equation of motion (88) is 



— — = tan sm Tj — tan sm n+i 

dt 2 2 + 



(92) 
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dt 
+ tan 



cos n ( cot Ki tan 



H+l 



esc cos Tj+i — tan 



— cot tan 



2 / 



K i-2 



CSC COSTj_i 



(93) 



From (93) we immediately compute 



dt 



5> = ° 



This is the discrete version of the conservation of total torsion H_ 2 in (22). We also note 
that in the continuum limit, (91) becomes (23). 



We now show that the equations of motion (|92|), (|93|) coincide with those of the LNS2 

(94) 



model. For this we introduce the following discrete variable, 



ipi = tan y e t " di 



where we have the (anti) symmetrized combination 



/i N 

5 2>-E 



Tk 



,fc=l 



k=i+l 



Note that (94) does not remain invariant under (82), (83). It should not be interpreted as the 
proper discretized version of the frame independent Hasimoto variable (16). Consequently 
we proceed by explicitly assuming the Frenet framing. 



A direct computation, using (90), gives 

dipi 



dt 

On the other hand, 

(l + M 2 )(Vm + ^ 



2 ii9 dK i 

sec — e ' — 
2 dt 



ztan— e l1?i — - 
2 dt 



sec 



2^ ^an^+V^+tan^V 



2 V 2 

By absorbing a factor of 2 in the definition of t and redefining 

ipi -»■ Vie" 

we obtain the LNS 2 equation 



This is a Hamiltonian equation of motion, with 

H = - ^ + = ~ 2 Y1 tan IT tan cos Ti +! 



and Poisson brackets 



= *(i + |^| 2 ) ^, 



{^,^-} = {^,^} = o 

These relations define the LNS2 model. Thus we have established the equivalence between 
the LHM and LNS2 models, by a direct computation in terms of the discrete Frenet frames. 



22 



X. DISCRETIZED RECURSION RELATIONS OF THE POSITIVE ORDER 



As an energy function for a discrete chain, the LHM (85), (86) engages only the tangent 
vector tj ~ Sj and consequently it is, as such, manifestly frame independent. But the 
combined variable (94) that relates LHM and LNS2 does not transform covariantly under 



frame rotations. Consequently it can not be interpreted as a discrete frame invariant version 
of the Hasimoto variable. Our goal, in the following, is to derive a proper i.e. manifestly 



frame independent discretized version of the continuum NLSE hierarchy (53) of conserved 



charges. These can be utilized as energy functions, in the case of piecewise linear polygonal 
chains. 

As the guiding principle we shall utilize solely the concept of invariance under frame 
rotations, in terms of the formal identification as gauge transformations: We demand that 
the discretized versions of the conserved charges should be manifestly invariant under the 
local frame rotations (80) i.e. they should derive from a gauge principle. They should 



go over to their continuum NLSE hierarchy counterparts, including the perfect derivative 
contributions, and smoothly in the limit where the bond length a between the neighboring 



vertices (79) vanishes 



We note that the primary emphasis here is on ensuring frame independence rather than on 
preserving integrability. The former is a necessary requirement for constructing meaningful 
energy functions for discrete chains, while the latter is a bonus. We presume that if the 
discretized recursion relations are carefully constructed so that frame invariance is preserved, 
these relations will also lead to integrability, if indeed the latter can be consolidated with 
the former. For this, we shall confirm that in a carefully executed continuum limit we 
indeed obtain the conserved charges of the NLSE hierarchy, including the perfect derivative 
contributions. In combination with frame independence, this is a strong indication that the 
discrete hierarchies we present are integrable. 



We introduce the following discrete version of the Hasimoto variable (15) 



(95) 



k=l 



Note that (95) is different from (94). The latter is not manifestly frame independent, while 
(95) is constructed to be so. We continue to identify the Ki = Ki+n as the bond angles (75). 



But from the conceptual point of view, there is now an advantage to consider the bond 
angles Ki as variables that have support at the corresponding i th vertex. 

Similarly, the variables Tk+i^k are related to the chain geometry, as torsion angles according 
to (76). Unlike the bond angles that reside at the vertices, the Tk+ik are now interpreted as 



variables that have support on the link that connects the vertices k and k + 1. 

As geometric bond and torsion angles, (k^, r i+1 ti ) relate to the discrete Frenet framing 
of the chain. Again, in the general frames, we extend Ki into a complex variable. The site 
dependent frame rotation that sends 



K i+ i ti = Ki ->> Kiexp{iA i+1 } 

Ti+li — > — H^i+l — 



(96) 



leaves the Hasimoto variable (95) intact (with Ai = 0). Note that in relation to (83), we 



have re-defined the torsion angles by scaling them with the bond length a. 
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We shall now introduce the discrete analogs of the densities u n (s) in (48). We denote 
them itf n [z]. We propose to evaluate the w n [i] from the following discretization of the (general 
frame) recursion relation (53), 



{w n [i 



W. 



,{i - 1] } + \Ki^2wi[i]w n ^i[i 



l<n 



The recursion starts with 



w [i\ 
wi[i] 



(97) 



(98) 



We have here carefully constructed (97) using the standard lattice gauge theory [42J dis- 
cretization of covariant derivative. We note that there are variants of this discretization [42J. 
In particular we may introduce 



{ w n [i + l}e iaT ^+ 1 - w n [i] } + XRiY^wi\i]w n ^[i 



l<n 



But for the present purposes it is sufficient to consider the version (97). We note that with 
the initialization (98), the quantities w n [i] transform according to 



iA,, 



l W n \l\ 



under discrete frame rotations. In particular, (97) is covariant under frame rotations. 
In analogy with (50), we introduce the following frame invariant quantities, 

N n [i] = RiW n [i] 



(99) 



and in the following we proceed by using the Frenet framing where We verify 

that in the continuum limit a — > the ensuing quantities (99) reproduce the densities (56), 
including the perfect derivative contributions. In combination with discrete frame rotation 
covariance, this will support the conjecture, that (99) are the proper discrete versions of 
the continuum conserved NLSE charges. We confirm this explicitly, for the four densities in 
([55])- 

In the case of Wi[i], it is clear that in the continuum limit the quantity N% [i] is reduced 



to 



lim Ni \i] 

n->0 



K 

T' 



(100) 



which coincides with the density Hi in (55). 
For w 2 [i], we get 



w 2 \t\ 



[wtli] - 



HI 



4 



In the continuum limit, this gives 



1~ tv 

lim N 2 \i] = 



i d 
8ds 



(« 2 ) 



(101) 



(102) 
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which is the conserved density ^(s) of NLSE hierarchy, including the perfect derivative 
contribution. 

The u> 3 [i] reads 

w 3 [i] = - { w 2 [i) - e~' MTl+1 ' l w 2 [i - 1] } + \KiWi[i)wi[i) 
ct 



a I a L 4 



4 J 



-tOTi+i . 



4 " 4 



IT 



In the continuum limit, 

lim iV3 [i] 



a->0 



1 // * / * / 2 , 1 2 2 i ^ 4 
■-KK TKK H T K H K 

4 2 4 4 16 



(103) 



/ 2 \ ' 

KK TK \ IK 



^ T 2 K 2 A 



4 4 16 



(104) 



This coincides with the density of the nonlinear Schrodinger Hamiltonian, together with the 



perfect derivative term in (55). 
For the fourth quantity, 



2 

w A [{\ = - { W3 [i] - e ~ WTi+1 ' i w 3 [i - 1] } + 2K i w 2 [i}w 1 [i] 



(105) 



In the continuum limit, we have verified that the ensuing charge N^i] reproduces rl±{s) in 
(55), inclusive the perfect derivative term. 



XI. DISCRETIZED RECURSION RELATIONS OF THE NEGATIVE ORDER 



We proceed to the discrete version of the negative order hierarchy (59). We propose the 
following manifestly frame rotation invariant discretization, 



(n+l) 



clk 2 



{ K 2 w„ n [i] - K 2 _ x e- iaT ^w. n [i - 1] } + A«?£«;_,[z]i«_ (n _ [i] (106) 

l<n 



As before, 



w [i] = 



and we initialize the recursion by setting 



w-i W 



This ensures that u>_( n+ i)[i] transforms covarianty under discrete frame rotations, 
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The higher order negative order quantities iV_ n [i] are defined as 



N. 



(107) 



In the following we proceed using the Frenet frames; by construction, the final results are 
frame independent. 

We have explicitly confirmed, that in the continuum limit this produces the correct Frenet 
frame densities in (61), including the perfect derivative contributions. For example, we 

iV_i = 1, 



immediately confirm that 

Similarly, for W- 2 [i] we have 
w- 2 [i] = 



CLK; 



w-iW 



<_ie" ,flT,+1 ''W-i[j 



1]} 



In the continuum limit, (107) gives immediately 



lim N_2 [i] = t — z(ln k)' 

a—>0 



(108) 



(109) 



which coincides with the density T-L- 2 in (61), including the perfect derivative. Similarly, 
we find that iV_3[i] produces the Proca mass term H_^{s) in (61) including the perfect 
derivative term: We have 



W-3 i 



2«J 



W- 2 \l\ 



i - 1]) + AK-w_i[i]w_i[i] 



We have also confirmed by explicit computation, that AL 4 [i] produces 1-L-^s) in (61 ), includ- 
ing the perfect derivative. Consequently, we conjecture that this is the case to all negative 
order terms. 

Finally, we note that if we introduce the discrete version of the Weyl transformation 



AL„[z] ^e^i\L n [z] = N e _ n \i] 



with 



6[i] = ln\Ki\ 2 

we find that in the continuum limit we recover the densities with positive order, 

a->0 



AL n [i] ^ Ni n \i] ~iV n [i] 
We have confirmed this explicitly, to order n = 4. 



H n (s) 



XII. ENERGY FUNCTIONS FOR DISCRETE CHAINS 

In [2S]"En]) the following energy function has been introduced and applied to model 
folded poteins as space filling piecewise linear polygonal chains, 

N N 

E = J2{- 2a i+i (J i + ba i J i+ c ( (X i- m2 )' 2 } + J2i dJi + qa i Ji + eJ i} (^v+i = 0) 
1=1 1=1 

(110) 
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The functional form of this energy function is justified by a naive discretization of the fol- 
lowing linear combination of charges H a (a = —3, —2, 1, 2, 3), that appear in the continuum 
NLSE hierarchy and its Weyl dual hierarchy; see equations (66)-(71) 



E 



ds 



(I) - w 



Aa 4 + /iVj + J ds{dJ + qa 2 J + eJ 2 } 



(111) 



Here the frame independent i.e. gauge invariant supercurrent variables are utilized. The 



energy function (111 ) has a natural gauge theory interpretation: The first line is the Abelian 



Higgs model in the supercurrent variables, with spontaneously broken potential term when 
fi 2 < 0. The second line displays the Chern-Simons term, the canonical momentum and the 
Proca mass term, respectively. When a naive discretization of the derivative contribution is 
implemented, 

da s 



ds 



(J 



2a- - 2a 



i+l 



we clearly arrive at (110). This discretization of the energy function (111) appears to 



correctly describe various properties of folded proteins 



Unfortunately, (110) also leads to a conceptual predicament. As the bond angle, the 



variable at in ( 110 ) is by construction a non-negative quantity, and in particular the first term 
in (110) does not remain invariant under a local transformation that send a^ — > —a^ at one 



given site. On the other hand, (111) is even in a and consequently in this continuum energy 
function, a(s) can be extended into a real valued variable. In the discrete case, we have 



previously used the Z 2 symmetry transformation (84) to extend the bond angle variable a. 



Ki from positive values into both positive and negative values. This symmetry transformation 



leaves the vectors tj intact, and as a consequence the discrete chain constructed from (78) is 



also intact. Consequently, from the point of view of discrete chains, it might be desirable to 
have a discrete energy function that remains invariant under this symmetry transformation. 



Unfortunately, (110) does not have this property. 



We note that the continuum version of the Z 2 transformation corresponds to a discontin- 



uous frame flipping by n at an inflection point of the curve; the continuum version of (84) 
engages both the unit step and the Dirac 5-function 



Our goal is to minimally improve (110) into a discrete energy function which is con- 



structed from the densities in the discrete NLSE hierarchy, so that the ensuing energy has 



the local Z 2 symmetry (84). In parallel with (111), we are instructed to try and proceed 
by using the discretized versions of the charges H a (a = —3, —2, 1, 2, 3). These charges have 
been derived in the previous two Sections, these are the quantities N a [i\; a G Z. These 
quantities are by construction fully frame independent and consequently invariant under the 



local symmetry transformation (84). Consequently we start by introducing the Z 2 invariant 
discrete combination 



a 3 N 3 \i) + a 2 N 2 [i} + aiA^i] + a_ 2 iV_ 2 [i] + ar_ 3 jV_ 3 [i] 



(112) 



By construction, this overlays the putative energy function. In general the charges N a [i] 
have both a real part and an imaginary part, 

N a \i] = ReN a [i\+iImN a [i\ 

and these are both independently invariant under the Z 2 transformation; we have observed 



(55), (61) that in the continuum limit, the imaginary parts are perfect derivatives, except 
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for a term in "H_3. As a consequence, whenever N a [i] is complex, we can interpret both the 
real part and the imaginary part, separately, as contributions to the energy. Accordingly, 
we account for them independently by substituting 



N a \i] -> ^ReNaii] + (3 2 ImN a [i\ 



P 



1,2 



(113) 



in the expansion (112). The ensuing real-valued quantity can be used as an energy function: 



According to the results in the previous Sections, its continuum limit contains the density 



in (111), in addition of terms that are perfect derivatives. In this way we can construct a 



Z2 transformation (84) 



discrete energy function that extends (110) in a manner which ensures invariance under the 



Unfortunately, (112) is quite elaborate in comparison with (110). But we have found that 



a discrete energy density which extends (110) in a manner which is invariant under (84) and 



yields the density in (111) in the continuum, can also be obtained by considering only the 



three quantities N a [i] with a = —2, 1, 2. 

We start with Ni[i]. It gives the following mass term contribution, we now work in the 
general frame where the complex 



(114) 



with m a parameter. 



We proceed to N 2 [i] in (101). We obtain the contribution 

JV 2 [i] ~ RiKi - /t i e" iaTl+1 ' ! K i _i 



(115) 



The first term can be combined with (114), it amounts to a redefinition of the a priori free 
parameter m. 

Following Section VII, we conceptually identify the second term as the standard U(l) 
gauge invariant kinetic term that is commonly introduced in lattice Abelian Higgs model 
[42J. For this, we simply interpret 

Kj ~ 4>i+l 

as the complex lattice Higgs field, and 

T i+l,i ~ Ai + ij 

the ensuing gauge field so that the second term is the familiar 



We write the second term in (115) as follows 



i?.p~* aT *+i>»K-. , 



exp{-iar i+1) i + In } 



ft , 



(JiOi-x exp{-iaJ i+ i )i } 



(116) 



where, in analogy with (65), we have introduced the gauge i.e. frame independent variables 



J i+1>i = r i+ i,i - i{arg(Ki_i) - arg^)} 
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The quantity (116) has both a real part and an imaginary part, that are independently frame 
rotation invariant. Consequently, in line with (113), we obtain two independent quantities 

(117) 



that we utilize to introduce the following contribution to discrete energy density 
E 2 [i] = y^{-2crjcr j _ 1 cos(a Jj+i.i) + sin(aJ i+M )} 



where the parameters are chosen for easy comparison with (110). 

Following Section VIII we proceed and extend the variables <7j from their a priori positive 
to arbitrary real values. For this we introduce the local Z 2 symmetry (84), that acts on the 
present variables at a given site k as follows: 



a 



J, 



— )■ 



for all j > k 



k+l,k 



+l,k 



(118) 



Unlike (110), the energy (117) remains intact under this symmetry, thus we can utilize it to 
extend a, in (117) from positive values to arbitrary real values. 
When we expand (117) in the bond length parameter a, we get 

« -loiOi-x + qOiOi-xJi+x.i + d^OiOi-x J? +li + 0(a 3 ) 



which is reminiscent of the first, fifth and second terms in (110), respectively. The difference 
is a higher order correction in the bond length parameter a. 



We now proceed to the negative order term iV_2[i]. From (108) we extract the following 
discrete energy density contribution. As in (113) we account for the real and imaginary part 
separately, and write 



2e o-i d Oi 

cos{aJ i+1 , i ) H srn(aJ; +M 



+ 



_ d 
a 

Qj - Oj-l 
a 



a a. 



(119) 



i-l 



sm(aJ i+ i t i 



2e 



cos(aJ i+ i 5 



dsm(aJ i+hi ) 



2e 



cos(aJ i+hi 



(120) 



Again, this remains intact under the Z 2 transformation (118) so that we are able to take <7j 
to be a real valued quantity. When we expand in the bond length parameter a, the first line 
gives modulo a constant 

w dJ i+1>i + eJf +M + 0{a 2 ) 



which is reminiscent of the fourth and sixth terms in ( 110 ) with corrections that are of higher 
order in a. The second line in (120) gives 

2e 1 



a 0-j_i 



+ 



Oj - Q-j-i 



(dJ i+hi + eJf +1>i ) + 0( 



where, in the continuum limit, the first term becomes a perfect derivative, 



1 Oi — o~i-\ d\na 
■ — — >■ 



Cj-l 



ds 
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and the second term vanishes in this limit, as 0(a) 



Finally, the quartic self-interaction of <7j arises naturally from (112), see (103), (104). 



Here, we propose to add it essentially ad hoc, as the square of the density (114); any even 



polynomial of is invariant under the discrete transformation (118) and as such admissible 
in energy. Most likely this is in violation of integrability. But in applications to proteins, 
it is desirable to have soliton solutions for which (discrete) translation invariance becomes 
broken, by a Peiers-Nabarro barrier |43j . 



version of the energy density (110) 



We combine (114), (117) and (119) to arrive at the following improved and Z2 invariant 



N N 

E = ^ I -2(jj(7j_ 1 cos (aJ i+ i ; i) + -a i a i - 1 sxn(aJ i+1; i)^ + J^. 
i=i a i=i 



2 2\2 

c (at - m L Y 



N 



i=l 



d a; 



a ovi 



sin(aJ i+1 , A 



2e Oi 
a 2 Ui-i 



cos(aJ i+u ) 



(121) 



It would be interesting to see, how (121) models folded proteins. 



XIII. SUMMARY 

In summary, we have addressed the problem how to utilize extrinsic geometry, to de- 
rive Hamiltonian energy functions for continuous and discrete strings that move in three 
dimensional space. The approach we have chosen, is based on a simple geometric principle: 
We insist, that the energy density must remain invariant under local frame rotations and, 
in the case of continuous strings, transform covariantly under reparametrizations. In this 
manner we arrive naturally to energy densities, that relate to the integrable hierarchy of the 
nonlinear Schrodinger equation. Furthermore, we have found that there is a Weyl dual to 
this hierarchy, that also relates to the energy densities that are relevant for strings in three 
space dimensions. We have argued that this additional hierarchy is also integrable, and we 
have confirmed the integrability explicitly in the case of the first few nontrivial quantities. 
In addition, we have established, by explicit computation, the equivalence of the two known 
integrable discrete nonlinear Schrodinger hierarchies. We have also shown, how a discretized 
generalization of the nonlinear Schrodinger equation that has been previously shown to de- 
scribe folded proteins, can be obtained as a particular limit of the quantities that appear in 
the discretized Zakharov-Shabat recursion relations that we have proposed. 

We hope that our approach can provide a systematic basis for the general description of 
both continuous and discrete string-like configurations in three space dimensions, in partic- 
ular space filling ones. 
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